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REPORT ON THE QUALITY OF THE 1955 COTTON CROP 


M. Lewin, S. HOROWITZ AND H. V. MusHAM 
Institute for Fibres and Forest Products Research, Ministry of Agriculture, Jerusalem 


ABSTRACT 


The 1955 Israeli commercial cotton crop consisting of Acala 4-42 was quite uniform from 
the spinners point of view. The fibres were of good quality, similar to the American fib- 
res of the same variety. The early picking period in the Northern part of the country yielded 
generally higher quality fibres with respect to the fibre strength, maturity and fibre length. 
The type of soil did not materially effect the fibre properties. All the results were submit- 
ted to the appropriate statistical treatments, which are briefly described in an appendix. 


INTRODUCTION 


The 1955 cotton crop was the first of commercial importance in Israel, and it was 
therefore decided to test a relatively large number of samples for the main fibre 
properties. 

Such tests are interesting for three main reasons: 

1. The assessment of possible differences in the properties of the fibres under 
different conditions might serve as a guide for the cotton grower in planning his 
future planting program. 

2. The comparison of the fibres grown under local conditions to international 
standards. 

3. The spinning mills require the highest possible quality and, perhaps even more, 
a uniform supply of cotton fibres. Differences in the quality of the cotton grown in 
various areas, on different types of soil and picked at different dates may be consider- 
able and present serious problems in production. Accurate information on the fibre 
properties may facilitate the choice of cotton from the proper lots and provide data 
on which the spinning mills can base an efficient program of blending of the various 
lots obtainable in the country. 

As only one variety of cotton had been planted — Acala 4—42 — in various 
regions of the country and on all types of soil, it is possible to assess the influence 
of these as well as that of the picking season on the properties of the cotton. The 
whole crop was grown on irrigated soil. 

The cotton was ginned on three gins, Lachish serving the southern, Param the 
central and Beit Shean the northern region, respectively. 
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The different types of soil were classiffied as light, medium and heavy. No cotton 
was grown on light soils in the Lachish region. It should be emphasized that this 
classification is a general one and does not yet take into consideration the soil 
properties as determined analytically. 

The first cotton was picked in August and the latest, of which samples were 
submitted to the Institute, was picked in January. Some cotton was picked also in 
February and March but no samples were submitted for testing. Most of the cotton 
was picked in September—October, a smaller part in November and only minor 
quantities from isolated fields were picked at later dates. Of the 183 samples tested, 
86 came from Beit Shan, 57 from Param and 40 from Lachish. 6 were picked in August, 
71 in September, 63 in October, 31 in November, 5 in December and 7 in January. 
35 samples were grown on light soils, 95 on medium soils and 53 on heavy soils. 


EXPERIMENTAL 


Of the total crop of about 10,000 bales, a 200 gram sample was drawn from each of 
183 bales and submitted for testing. The samples were drawn by the Cotton Marketing 
Board so as to represent proportionally the regions, types of soil and picking periods. 
The samples were cut out from the bales after being pressed at the gins in accordance 
with the standard practice demanded by the grader, were first judged by the grader 
and afterwards submitted to the Institute; here they were blended on a cotton 
blender, designed and built by the Institute. 


The strength of cotton fibres was determined by means of the Pressley Strength 
Tester in accordance with A.S.T.M. Designation D1145—53T.1 


For determining the length of cotton fibres the array method was used. The 
instrument employed for the tests was the Suter-Webb in accordance with ASTM 
D414—S4T. 

The fineness of cotton fibres was tested by means of the Micronaire in accordance 
with ASTM D1448—S54T. 


The maturity of the cotton fibres was determined by means of the Causticaire 
method as used by the United States Department of Agriculture. 


RESULTS AND DISCUSSION 


In Figures 1—5 the changes of the tensile strength, maturity, fineness, U.Q.L., grade 
and grader’s length with picking time for the three types of soil are summarized. 
Each point on these figures represents the calculated average of the values deter- 
mined experimentally on each of the samples of the appropriate group. 

In Figures 7—12 the changes in each of the properties investigated with the picking 
time for the 3 gins are summarized. 


In the following the individual properties are discussed with reference to the 
statistical significance of the results obtained. 
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Figure 1 
Changes in the tensile strength (Pressley Index, 
zero gauge) with the picking time of cotton fibres 
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Changes in grade (determined by Grader) with 
picking time of cotton fibres grown on various 
soils 
Legend as in figure 1 
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Figure 2 
Changes in maturity (Causticare method) with 
picking time of cotton fibres grown on various 
soils 
Legend as in figure 1 
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Figure 6 
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Figure 7 
Changes in Pressley Index (zero gauge) with 
picking time of cotton fibres from 3 different 
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Strength 
The average strength of the 183 samples gave a Pressley Index of 8.42 -+ 0.06. This 
compares with about 8.2 obtained for Acala 4—42 in the U.S.A. 

All of the factors studied — soil, gin and month of picking — showed some sig- 
nificant effects on the strength of the fibres, expressed as Pressley Index. 

Whatever the effect of the month of picking on the strength of the fibres was 
significant (Beit Shean, heavy and light soils) the P.I. decreased in the latter part 
of the season (see Table I). A similar trend was observed with the 1955 cotton crop 
in 5 out of 12 cotton growing states in the United States, while in none of the States 
a statistically significant trend of the opposite sense could be observed2. 


TABLE I 
Pressley Index by month of picking 


Aug. Sept. Oct. Nov. Dec. Jan. 


P. I. Beit Shan, Light Soil — 9.09 8.61 8.07 8.48 8.29 
P. I. Beit Shan, Heavy Soil 8.76 8.68 8.28 8.5 8.21 ~ 


In two other instances, Beit Shean and Lachish, medium soil, where a sufficient 
number of samples from the Israel crop were available to test for the significance 
of the effect of the month of picking, one series showed a decrease, the other an 
increase in the P.I., both being statistically non-significant. 

It appears (see Table II) that for each of the seasons of picking for which the effect 
of the gin on the P.I. could be studied the P.I. of cotton from the Beit Shean area was 
among the higher ones, independently of the types of soil. Thus, for instance, samples 
from Beit Shean showed P.I.’s of 9.09 and 8,63 (medium and heavy soil respectively 
against 8.51 (Param) and 8.14 (Lachish); in October 8.30 (both Beit Shean and 
Param) against 8.20 (Lachish); in November 8.02 (both Beit Shean and Lachish) 
against 7.96 (Param). 

Both in September and October, Param gave significantly higher P.I. (8.51 and 
8.30) than Lachish (8.14 and 8.20) but in November the position was reversed 
(Lachish 8.20 — Param 7.96) (Table II). 

Samples taken from light soil consistently showed relatively high P.I. No such 
conclusion can be drawn with respect to a difference between heavy and medium 
soils (see Table IT). 

Fibre Length (U.Q.L.) 
The average U.Q.L. for 183 samples was 1.089 +0.0052 or 13/32 app. This corresponds 
to U.S.A. results for Acala 4—42. 

The U.Q.L. decreased consistently as the season progressed. This trend was 
statistically highly significant in one instance (Beit Shean medium soil) and slightly 
significant in one other case (Beit Shean light soil). A decreasing trend in fibre length 
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from early through midseason to late crops was also observed with the 1955 crop 
in the U.S., but there the trend was less conspicous than it appears to be in Israel. 


No consistent effect of the gins on the U.Q.L. could be observed (see Table II): 
Beit Shean showed relatively high U.Q.L.’s in September and low ones in November, 
but in October low U.Q.L. on light and heavy soil but high U.Q.L. in medium 
soil. With respect to Param and Lachish differences in U.Q.L. are even less consistent, 
as they are always significantly affected by the type of soil too. However the effect 
of the soil itself is also not uniform. When heavy soils from Param yielded significantly | 
short fibres (U.Q.L. = 1.093) but at the same time heavy soils from Beit Shean 
yielded short fibres (U.Q.L. = 1.042). Light and medium soils could be found to. 
yield most inconsistently both the longer and shorter fibres, according to the region | 
in which they had been grown. 


Fineness | 

The average Micronaire Number for the 183 samples was 4.35 -++ 0.031 comparing | 
with 4.2 obtained in the U.S.A. 
The effect of the month of picking on the Micronaire Number could be tested | 
for only four combinations of gin and soil. In two cases the trend was upward and | 
in two cases downward. No conclusions could therefore be drawn on the effect of the | 
month of picking on the fineness. It might be mentioned that with the 1955 crop in | 
the U.S. early samples tended to be coarser and late samples to be finer?. | 
With respect to the effect of the gin and the soil on the Micronaire Number no | 
clear cut effects could be observed either. | 
The only statistically significant feature was an “interaction” effect between gins — 


and soils, i.e. that for each of the gins the effect of the soil was of a different nature. — 


Maturity 


It should be mentioned that as the data of sowing were not identified for the samples, 
there is not necessarily a correlation between the date of picking and the age of the 
plant. 

The average for all the samples gave a maturity of 80.0 + 0.8. This equals U.S.A. 
results. 

The Maturity decreased consistently as the season progressed, the trend being 
highly significant in two instances (Beit Shean heavy soil, Lachish medium soil). 
Similar results were observed in the U.S.! where in eight out of twelve cotton growing 
states fibres from early samples were slightly more mature than those from midseason 
and late. 

Beit Shean showed consistently higher maturity (September 81.9, November 80.28) 
than Param and Lachish (September 80.81 and November 76.76). No significant 
difference in maturity between samples from the latter two gins could be observed. 

The effect of soil on maturity was found to be non-significant. 
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Grade 


The grades were determined by the Marketing Board’s Grader, Mr. Gefen. 

To enable computations, numerical values were given to the various grades. The 
average for all the samples is just short of strict middling. 

The difference between the picking periods was found significant: August, Sep- 
tember and October yielded an average of better than strict middling while November, 
December and January gave an average of just better than strict low middling. 

The difference between the soils was found not significant. 


CONCLUSIONS 


The 1955 crop was quite uniform from the spinners point of view. This is quite 
logical as only one variety was grown. The cotton was of good quality similar to 
that of the American Acal 4—42. 

The early picking period yielded generally higher quality fibres with respect to 
P.I., Maturity and U.Q.L. It is important to note that the first rains fell in the middle 
of the picking season so that the earlier samples were picked before and the latter 
after the onset of the rains. Gin I gave generally better results than the other two 
gins, while no overall effect of difference of the soils could be observed. 

Fineness, expressed as Micronaire number, averaged 4.35. Each of the three 
factors showed statistically significant effect on the fineness but of inconsistent 
direction. 

The maturity index determined by the Causticaire method averaged 80.0. The more 
mature fibers were found in the earlier part of the season and in the Beit Shean area. 

The average grade as determined by the official grader was just short of strict 
middling. The earlier part of the season yielding a grade somewhat better than 
strict middling while the latter part of the season, just better than strict low middling. 


APPENDIX I 


STATISTICAL ANALYSIS 
SEZIACKS 


The statistical analysis of the data was designed to test the significance of the respective 
effects of gins, soils and data of picking on various characteristics of the fibres. 
Results were classified according to these three factors, and graphical representations 
of sample means of each variable against data of picking were plotted for each 
combination of gins and soils. 

As these diagrams showed certain trends, the significance of such trends against a 
null hypothesis of randomness was submitted to a nonparametric test. For this 
purpose a test based on Kendall’s rank correlation coefficient*+ was used and the 
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probability of obtaining, in the absence of any trend, a value exceeding the observed 
one was used as criterion. 

The results of this procedure are given in Table III. Seven out of sixteen series 
tested by this procedure revealed statistically significant trends. As the date of picking 
could thus be shown to be a significant factor, the effects of gins and soils upon 
the measured characteristics had to be tested for each month of picking independently. 


TABLE Ii 


Effect of season on fibre properties 


Fibre property Gin Type of soil Direction of trend Significance 


Pressley Index Beit Shean Light decreasing + 
Medium decréasing _ 
Heavy decreasing +++ 
Lachish Medium increasng —_ 
Micronaire Beit Shean Light decreasing — 
Number Medium increasing — 
Heavy decreasing +++ 
Lachish Medium increasing — 
Maturity Beit Shean Light decreasing _ 
Medium decreasing — 
Heavy decreasing +44 
Lachish Medium decreasing +44 
Upper Quartile Beit Shean Light decreasing + 
Length Heavy decreasing +++ 
Medium decreasing — 
Lachish Medium decreasing — 
Notes. +++ = significant at a level of P < 0.05 
++ = significant at a level of 0.05 < P < 0.10 
a = significant at a level of P = 0.10 


= non significant 


increasing means low values at early season high values at late season. 
decreasing means high values at early season low values at late season. 


The months of picking chosen for this analysis were, August, September, October 
and November, too few samples being available from the other months. 


In order to test whether the effects of gins and of soils are significant, an analysis 
of variance according to a two way classification was carried out according to a 
model which includes an interaction factor between the two factors. This model is: 
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Xie = + ae + By + yg + eng 
where a; stands for a gin factor at a level i; 6; for asoil factor at a level j, 4 for an 
interaction between a; and §; and ey, is a random error variable, normally distributed 
with zero mean and variance o2;;. The assumption about the distribution of ej, was 
confirmed by means of a nonparametric test concerning the Kolmogorov Smirnov 
d-statistic® which is independent of the parameters of the normal distribution under 
investigation, this test confirmed our assumption. 

The most serious problem which the analysis of variance had to encounter was 
aroused by the differences between the size of the various samples, which introduce 
a lack of an orthogonality of the factors a; and f;. Hence, a straightforward procedure 
of analysis of variance for a two way classification could not be applied and a more 
cumbersome procedure of fitting the constants a; and f; according to the least 
squares principle had to be followed. 

However there remained two circumstances where even this procedure of analysis 
could not be followed, namely: 

(a) when the sample estimates of the within class variances (02;;) were significantly 
different on the basis of Bartlett’s test (Bartlett’s test was in fact, significant at the 
5% level only in two instances, micronaire number and U.Q.L. in September. In 
all other cases it was non-significant event at the 10% level). 

(b) where many combinations of soils and gins were missing so that non-systematic 
cross classification could be obtained. In these cases, as well as whenever the analysis 
of variance had shown significant effects of at least one of the factors, two methods 
of separating sample means into homogeneous groups were. applied. 

These methods are due to Duncan} and Tuckey’. 

At last certain remarks should be made about the choice of the significance level 
of the tests. It is most common to take P= 0.05 as a significance level. However, 
for a given number of degrees of freedom, if Pis decreased the power of the tests is 
also decreased. Hence, in order to keep the tests sensitive it was decided that in 
cases where the number of degrees of freedom is small P would be increased to 0.10. 

The results of the statistical analyses are summarized in Table II. 
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SUMMARY 
Samples of Acala 4—42 cotton fibres grown under irrigation and representing the 
1955 commercial crop were drawn from the three main regions, three types of soil 
and throughout the whole picking season and tested for their basic properties. 
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The fibre strength expressed as Pressley Index averaged 8.42, the highest values 
occurring in the early part of the season in the Beit Shean area and on light soils. 

The fibre length (U.Q.L.) averaged 13/32”. The longer fibres were obtained in the 
early season, but no clear effect of the region or the type of soil on the fibre length 
could be observed. 
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REFLECTION AND TRANSMISSION OF INFRARED RAYS BY 
INDIAN TIMBERS 


D. NARAYANAMURTI AND U. SANJIVA 
Forest Research Institute, Dehra Dun., India 


The drying of veneers by means of infrared radiation was discussed in two earlier 
papers!,2. It was thought of interest to determine the transmission, reflection and 
absorption coefficients for infrared of various Indian timbers at different moisture 
contents and temperatures, as only meagre data on the subject is available in the 
literature. 

Prat? investigated the permeability of wood to infrared rays by inserting the wood 
between an infrared sensitive photographic plate and the source of radiation (150W 
lamp at 2 m distance), an exposure period of 6 hours being given. He found Larix, 
Abies sp., Picea sp., Alnus and Tilea sp. the most permeable of the lot, ie. they 
could transmit the rays through a thickness of 5—7 mm. Pinus sp., poplars and Acer 
sp. showed medium permeability (3—4 mm) while, Carpinus, Betula, Fagus, Ulmus, 
Robinia, Quercus, Prunus, Juglans and Pterocarpus were least transparent (0.5—2mm). 
Rawling4 and Clark5, using 3 mm specimens, also classified the timbers into four 
groups. 


TABLE I 
Very Transparent Moderately transparent Almost impermeable Impermeable 
Pinus cembra Pinus Tectona grandis Juglans regia 
Pinus Eucalyptus Swietenia Terminalia bialata 
Acer pseudoplatanus Quercus Diospyros 
Buxus 
Fagus 


Deribére® studied reflection coefficients at three ranges of moisture content and his 


results were as follows: 
TABLE II 


Reflection % of perfect reflection 
Moisture content range 


Poplar Beech Oak 
ee RW reks. TeRe Wialls. TRS Wale 
0-2 38.5 68 31.5 65 28.5 57 
10-15 30.5 56 29.0 54 28 47 
30-40 19 69 19 63 18.5 60 


The present note records the results of tests with some Indian species. 
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SPECIES OF WOOD INVESTIGATED 


40 species of Indian timbers varying in density from 0.2 to 1.3 g/cm} were investiga- 
ted in the air dry condition for the reflection coefficients, and several were studied 
in the wet state. The surfaces of the specimens were carefully prepared to the same 
degree of smoothness. 


METHOD ADOPTED 


For measuring the reflection, transmission, etc. a G.E.C. lamp (115 volt) with 
external parabolic reflector was used. The wood, mounted on a stand was at the 
focus of the lamp radiation. The reflected light was measured with the help of a 
Moll thermopile. For measuring the voltage of the thermopile, in the earlier experi- 
ments a milli-voltmeter was used. Later a Cambridge potentiometer with a galvano- 
meter with an accuracy of 0.01 mm was employed. For some experiments a Zernicke 
galvanometer was also used. Usually the rays impinged on the surface at an angle 
of 30°. In a few experiments the influence of the angle of incidence was also investi- 
gated. Measurements were taken with both rock salt and glass windows for the 
thermopile. For penetration studies the same set up was employed. However, the 
rays had normal incidence and the thermopile was placed behind the specimen. 1 cm 
thick specimens were usually taken and after each reading they were carefully thinned 
down till there was just some response in the thermopile. The influence of thickness 
penetration was also studied in a few cases. The variation of penetration with moist- 
ure content was investigated with veneers of some species at various stages of drying 
The voltage of the lamp was kept constant by means of a variac transformer and 
voltmeter. 


RESULTS AND DISCUSSION 


The percentage reflection of infrared radiation varies from 21.7% of the bright 
reflector in the case of ebony to 41.2% with Cryptomeria (see Table III). For the case 
of the timbers studied it appears to decrease with density (Figure 1). This is also 
indicated by the results of Deribére 
where poplar has the higher reflec- 
tion. The results of some tests in the 
wet condition are given in Table IV 
below. 

It is seen that the percentage of 
reflection decreases with increased 
moisture content. This is to be ex- 
pected as water is a good absorber 
of infrared radiation. The results are 
o2 04 6G O8 TO 1? similar to those of Deribére but in 

pay ae the case of Tetrameles nudiflora the 
Figure 1 effect is minimum. 
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TABLE III 
Relative reflection coefficients of Indian woods for infrared 


Sp) et. 74 Of So. gt.) Sor 
bright bright 
reflector reflector 
Artocarpus integrifolia (jack) 0.64 34.9 Prunus paddam 0.60 36.0 
Artocarpus hirsuta (aini) 0.64 36.4 Pentace burmanica (thitoha) 0.80 31.7 
Abies sp. (fir). 0.44 37.0 Pterocarpus macrocarpus 
Albizzia procera (siris) 0.64 39.4 (padauk) 0.88 32:2 
Cullenia excelsa 0.59 36.2 Pinus longifolia (chir) 0.44 38.1 
Cupressus torulosa (cypress) 0.57 36.6 Picea morinda (spruce) 0.57 39.8 
Cedrela toona (toon) 0.67 34.0 Shorea robusta (sal) 0.91 212, 
Cedrus deodara (deodar) 0.54 35.0 Schima wallichii (chilanuni) 0.74 33.8 
Calophyllum sp. (poon) 0.68 35.8 Shorea assamica 0.62 36.6 
Canarium sp. (dhup) 0.57 40.7 Swintonia floribunda 
Carapa moluccensis (kyana) 0.99 32.0 (taungthayet) 0.55 39.2 
Cryptomeria japonica (suji) 0.29 41.2 Swietenia macrophylla 
Diospyros sp. (ebony) 1.31 PAI (mahogany) 0.52 39.9 
Dalbergia cultrata (yindaik) 1.04 DTI Tamarindus indica (tamarind) 1.00 Qe 
Dalbergia latifolia (rosewood) 0.78 30.1 Terminalia paniculata (kindal) 0.86 31.0 
Dalbergia sissoo (sissoo) 0.82 34.4 Terminalia bialata (chuglam) 0.57 37.6 
Dichopsis elliptica (Pali) 0.77 33.8 Terminalia tomentosa (laurel) 0.97 28.9 
Hardwickia pinnata 0.71 B51 Taxus baccata (yew) 32.4 
Holoptelea integrifolia (kanju) 0.65 S73 Tectona grandis (teak) 0.67 34.8 
Juglans regia (walnut) 0.53 34.6 Terminalia myriocarpa(hollock)0.52 38.8 
Lagerstroemia flosreginae 0.76 35.6 Tetrameles nudiflora 0.40 39.9 
Mangifera indica (mango) 0.46 36.8 Xylia dolabrifomis (pyinkado) 1.09 26.6 
Michelia champaca (champak) 0.56 38.1 Zanthoxylum rhetsa 0.43 38.8 
Millettia pendula (thinwin) 1.04 26.0 
TABLE IV 
Species % reflection 
Dry Wet 
Tetrameles nudiflora 34.91 19.42 
Cryptomeria japonica 39.20 8.17 
Dalberga sissoo 40.57 21.88 
Eleocarpus sp. — 17.92 
Dysoxylum malabaricum — 22.92 
TABLE V 
Species Maximum thickness Species Maximum thickness 
(mm) (mm) 
Cryptomeria japonica (suji) 7.6 (9.2)* Homalium sp. 3.4 
Pinus longifolia (chir) 5.6 Terminalia myriocarpa 
Picea morinda (spruce) 5.6 (hollock) 3.0 
Cedrus deodara (deodar) 6.0 (8.5)* | Gmelina arborea (yemana) 3.8 
Taxus baccata (yew) R 4.3 Zanthoxylum rhetsa Bi 
Taxus baccata (yew) T Del Dipterocarpus sp. (gurjan) 4.1 
Calophylium sp. (poon) 2.0 Ochroma sp. (balsa) 7.6 
Terminalia bialata (chuglam) 2.0 (4.5)* | Mangifera indica (mango) —2.0 
Canarium sp. (white dhup) 3.0 (4.5)*  Tectona grandis (teak) — (6.2)* 
Artocarpus integrifolia (jack) 3.5 Dalbergia sissoo (sissoo) — (4.0)* 
Artocarpus hirsuta (aini) 212, Juglans regia (walnut) — (4.7)* 


* Readings taken with a Zernicke galvanometer. 
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A few experiments were carried out on the effect of the angle of incidence on reflection. 
Angles of incidence of 15, 20, 25 and 35° were employed. In general, the percentage 
reflected seems to be somewhat less at 35°. This is being further investigated. The 
influence of the nature of the surface is also being studied. 
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Figure 2 
Zanthoxylum rhetsa 1/32” 
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Figure 3 
Mango (air-dry) 


The results of penetration studies are 
given in Table V above. 

The above results indicate that generally 
conifers are more easily penetrated than 
hardwoods. The results of Prat} and 
Rawling4 and ClarkS indicate teak, wal- 
nut and chuglam to be rather imperme- 
able. Our results also indicate chuglam to 
be least permeable. Among the hard- 
woods balsa seems to be very permeable. 

The variation of permeability with 
moisture content and thickness is indi- 
cated in Figures 2, 3 and 4. In a few 
cases a slight increase is noted again at 
low moisture contents. 
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Figure 4 
Balsa (air-dry) 
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THE COMPLETE STRESS-STRAIN LAW IN INFINITESIMAL ELASTICITY 


M. REINER AND L. RINTEL 
Technion—Israel Institute of Technology, Haifa 


ABSRTACT 


It is shown that in infinitesimal strain, second order effects will be present when there is a 
rotation of the principal axes. The complete description then requires two new moduli in 
addition to the two of the classical theory. 


1. The elastic state of a solid body is described by the tensor of strain (e), to which 
the tensor of stress (s) is related in a single-valued function 


s = F(e) (1.1) 
or vice versa 


e = Xs) (132). 


In the classical theory of isotropic material, the first equation is postulated in the 


linear form 
Sim = Ar, Ey Om + 2 Elm (1.3), 


where €1m is the tensor of “infinitesimal” strain, 2, and mw are constant and éy is the 
infinitesimal cubical dilatation. 

2. The concept “‘stress” is unequivocally defined through the primitive concepts 
of “force” and ‘“‘area-vector’’. This is not the case with regard to the concept “strain”. 
To define the latter we have recourse to the primitive concept “displacement gra- 
dient” (I"). Let 1x be the “‘initial”, x; the “final’’ position of a particle. It is assumed 
that the coordinates x; are singlevalued and continuously differentiable functions 


of the 1x, and vice versa or 


xt = f(x) 5 1x = FOr) (2.1). 
The displacement is then defined by either of the two expressions 
uy = x1 — F(x1) 5 U = f(x) — ix (2.2); 


and the displacement gradient by 


Tea == Oui/OXm = O1m = OF (x1)/Oxm (2.3) 
tml’ = Ou/Omx = Of(x)/Omx — 1md ; 


where Om is ‘‘Kronecker’s delta’’. 
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constitutes a “measure of strain”. However, this measure is not unambiguo 
because any isotropic symmetrical tensor function of J’, such that it vanishes when 
represents a rigid displacement, is a suitable measure of strain. 

It should be noted that the displacement gradient is a dimensionless number 
“infinitesimal” in aaoaesnel elasticity” refers to ihe strain. - strain ed infiinii 


lations carried out in a problem, be neglected when added to unity. For instance | 
cylinder of length /, if made of steel, can be strained to a length ] = Ip (1 +0.001)) 
so that (J—/9)//g = 0.001, and not more. Such strain can be considered as infinid 
tesimal. 

3. The classical theory of elasticity is “infinitesimal” in the sense mentioned be 
fore. However, its procedure “begs the question’. In classical theory of elasticity, 


the infinitesimal strain « is defined by | 
Elm = 4(yim + Vm) G1] 


As will be shown in the present paper, the procedure of the classical theory thud 
suppresses second order terms which may arise even when the displacement gradient 
is infinintesimal, because infinitesimal second order terms can be neglected only 
when the first order term is present, and in comparison with it. For instance Love 
in his classical treatise defines finite strain by a measure which has as its components 
Cxx = Ouz/Ox + 4[(@uz/Ox)? + (Cuy/dx)2 + (Guz/0x)? | etc., and proceeds to argue 
that “the quantities ez... which were called ‘components of strain‘ in previous 
articles are sufficiently exact equivalents of ezz, when the squares and products of 
such quantities as u/0x are neglected’. However, consider the case when uz = co.) 
Then the first term of Love‘s complete expression is reduced to zero and the term: 
in brackets to (@uy/0x)* + (éu;/0x)2. This, as has been emphasized, by one of us 
cannot be neglected as “infinitely small of the second order“ because the first orders 
term is absent. “Infinitely small of second order‘‘ when there is no “‘first order“, 
is meaningless. The quantity ezz = 4[(@uy/@x)? + (éu./dx)2], if not vanishing,, 
cannot be neglected against zero in comparison with which it is infinitely large. 
There is, of course, nohing to prevent aninvestigator to disregard such a term provided 
he knows of its existence. If he does so without examining the consequences, such 
procedure may be dangerous. He can be likened to one driving a car on the street 
of a Dutch town at the side of a canal. It may be of no importance whether his forward 
velocity is 50 or 55 miles per hour, and he may safely neglect the additional 5 miles, 
but if his velocity component in the sidewise direction is of that “negligible“ order, 
this will land him into the water. 

The correct procedure is to solve the problem for finite strain assuming finite I”, 
and to carry out the limiting process with admissible approximations assuming 
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infinitesimal y in the solution. The result will reveal second order effects entirely 
masked by the classical procedure of linearizing at the head, effects which are ame- 
nable to experimental verification and have actually been observed. As shown in 
the example given before, second order effects may appear when the fiirst order 
of classical infinitesimal elasticity is absent. The method for discovering such second 
order effects will therefore consist in solving a problem in accordance with Eqs. 
(3.1) and (1.3), and looking for points and directions in which some em vanishes. 
Solving the finite problem for the same point and direction, and going over to infi- 
nitesimal displacement gradients may then reveal second order effects. 

4.A number of measures of finite strain have been postulated so far. When we 
consider the strain - ellipsoid with half-axes A(t), where ¢ is a running index for 
i,j,k, any function of A(t) which vanishes for A(t) = 1, when there is no deformation, 
will form a suitable measure. In terms of A(x) the lollowing basic measures have 
been postulated by Cauchy, Green, Hencky, Almansi and Swainger in this order.* 


: 
C=1-1 
G — 4022 — 1) 

t (4.1) 


H — Ina = — In(1/A) 
4 = 4(1-1/22) 
S=1-1/A 
Other measures result from combinations of these. For instance a measure proposed 
by Wall4 
Mi = 4a—1/) = 465) (4.2) 
and another measure proposed by Wall5 
W = (42—1/2) = A, Si) (4.3). 
The half axes J in terms of the different measures are shown in the first column 
of Table I. Let i,j,k be the principal axes, then the ratio Vo/V, where Vp is the initial 


volume of a portion of the body under consideration, and V is the volume occupied 
by the same particles in the strained position, is given by 


V/[Vo = Ai Aj Ak (4.4) 


with ey, the cubical dilatation, 
ey = (V—V)/Vo = V/IVo — 1 (4.5). 


* For references compare Truesdell. 
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This ratio for the different measures is shown in the second column of Table I, 
with I, IJ, III indicating the first, second and third invariant of the strain tensor, 
respectively. 


TABLE I 
Half axes of the strain-ellipsoid and volume-ratios for different strain measures 


1 2 
Measure A ev 
Cc 1l+e I++ 1 
G (1 + 2e)'/2 (1 + 21 + 411 + 8iip'/2—1 
H exp(e) exp(I) — 1 
A (1 — 2e)—I/2 (i — 21 + 411 — 8i)—'/2—1 
S (i — e)— dag—-1+0-—-)—!—-1 


For a general theory of elasticity we would need tensor expressions for these 
different measures. They have been derived by Murnaghan? for G(e) and A(e). Hanin 
and Reiner? have shown how to derive such expressions for any other measure, 
and more especially for H(e). However, this is not imperative for our argumentation. 
What we need here are expressions for the two cases of 


A) pure strain 
B) simple shear, 


These expressions can be derived by elementary methods following Love}. 


5. A) Pure Strain is one in which the principal axes do not rotate. In this case 
they can be chosen as axes of reference. We then have 


Ai => Ox; | 04x = 4X4 
1/Ag = O4x/0x4 == aX, } OD 


and from (2.3), introduing for indices / and m the same index i, which makes 6; = 1 


A= Tr 
y= Ver } G2 


Eqs. (4.1) now bring forth values entered in the first line of Table II. 
Now let J’ be infinitesimal, which we indicate be writing for it y. Then from 


L+iv=1/l—“) =1lt+nwty? +yb4+... (5.3) 
there follows 
ey (5.4)* 
and all measures are reduced to the classical measure 


= Yi (5.535 


* Note that this relation is valid only when no second or higher order terms are present. 
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TABLE II 
Principal pure finite strain 


Cc G 
Se oe | eee ee Tae ae eee ee Ee ee Pre Pee ee eee a 
1 e: it Has il 
Zeer a—1 4 (a2 — 1) 
3 e; b-1 4 (62 — 1) 
4 ex c—1 $2.24) 
Saale a+b+c—3 4 (@ + b2 + c2 — 3) 
6 Il (a—1)(6—1)+(6—1)(c—1)+ 4 [@—1)@2—1)+ 
a-Kc—)iia—1) +(b2—1)(c2—1)+(c2—1)(a2—1)] 

i bil (a—1)(b—1)(c—1) © 1/3(a2— 1)(b2—1)(c2—1) 

A, A, 
ares Ini+:D) — Ind — I) 
Denes Ina —TIn 1/a 
Bites In b —TIn 1/b 
4 e& Inc — In i/c 
a IE In (abc) — In (1/abc) 
6 Il (In a)(In b)+(in b) - (In 1/a)Un 1/6) +n 1/6)(n 1/0) + 

+ (Inc)+(n cin a) +(In 1/c)(In 1/a) 

7 Ile (In a)(In b)(In ©) —(In 1/a)(In 1/b)(In 1/c) 

A S 
1 e L:(1— L%/2) rT; 
2 & 4 (1 — 1/a?) 1 — l/a 
acs 4 1 — 1/62) 1 — 1/b 
4 ex 41 — 1/c2 1—I1/c 
Suh 4 (3 — I/a2 — 1/b2 — 1/c?) 3 — 1/a — 1/b — 1/c 
6 Ji. 4 [G1 —1/a2)(1—1/62) +1 — 1/6?) - (1—1/a)0 —1/6)+( —1/6)0 —1/¢e)+ 

+ (1—1/c2)+ (1 —1/c2)(1 —1/a?)] +(1—1/e)1—1/a) 

7 Ide 1/3(1 —1/a2)(1 —1/b2)(1 — 1/c?) (i—1/a) —1/6)0 —1/c) 


For instance, homogeneous pure strain is given by 


Kee Gh) ak j= DGX). GX Ea) (5.6) 
(eo =a - Cayeyx = 1/6. XP) iar = We-x) 


where the a, b, c are constants, so that from (2.3) 
W' =a—1 yak sb) 5 el == c—] (5.7) 
Ty =1—Il/a ; FF; =1—1/b ; 1, = 1-1 /e ; 


and from (4.2) 
fp ee AR EG = (5.8). 
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The principal strains are shown in lines 2, 3 and 4 of Table II, followed by the 
invariants of the tensor. Going over to infintesimal strain, when I » y, we find 
for instance for &, introducing from (5.2) for lim a =1 + ya, that $(a2—1) = 
va(1+vyal2) & ya, when we neglect y against unity. In this manner all measures are 
reduced to the same expressions, and we find the values entered in Table III. 


TABLE It 
Principal infinitesimal pure strains 


& &j Ek I . II é III i é, 
Ya Yo Ye YotYot+Ye Payot povet+YcVa Ya Yo Ve YatypotYe 


The principal strains are accordingly of first order only, and are the same as in the 
classical theory. 

We now consider two cases of homogeneous pure strain, namely (a) cubical 
dilatation, (b) linear dilatation. Their values can be found from those shown in 
Table III by specialisation. 


6. (a) Cubical dilatation is given by 
x VaGx) , x = Waesws (6.1) | 


where / stands for any index i,j,k or x,y,z. We accordingly find the values shown | 
in Table IV. | 


TABLE IV 
Strain-components in infinitesimal cubical dilatation 


(b) For /inear dilatation in the direction k we introduce the longitudinal extension 
yi for yx, and the transversal contractions — y; for ya and yp, and find from Table III 
the values entered in Table V. 


TABLE V 
Principal strain-components in infinitesimal linear dilatation 


& &3 Ek qe I, il, &, 
—yt —yt yu yi—2yt —yi(2yi— yt) yr yi yi—2yt 


7. We now go over to (B) simple shear as treated for instance by Love! (Art.37). 
This is kinematically given by 


X1 = (1x) +L (x2) , x2 = 2x , x3 = 3x (7.1) 
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from which 


Py2=pF=L (7.2) 


while all other components vanish. Introducing 


r=2 tan a (7.3) 
Love calculates 
A; = (1+sin a) / cosa, Ay = (1—sin a) / cosa, Ay = 1 (7.4). 
These equations yield the expressions shown in the following Table VI. 


TABLE VI 
Principal Strains in Simple Shear 


Cc G H A S 


1 e& 4 (/+2)—1 P/4-(WV4+Tl) mtk(V+D)]) L/4-(V-l) -4(V/-TP)+1 

26 $(V-P)-1 -T4-(WV-T) mV -P)1 -T/4-(V4P) -$(V4F)41 

BREIL? SoKse2 r22 0 27/9 Daas 

4. 2-vV =i] n[t(/—-D)]- —L2/4 2-/ 
‘n{h(Vt+l)] 

5 J =VI2+44 er=0 II.=0 &o=0 


However we want to know the components of the strain-tensor in respect to the 
system x, y, z in the final position. For this we make use of a relation proved by Love 
at the same place, namely that the directions of the principal axes i and / are the 
bisectors of the angle 7/2 + a with the x— axis, and that the angle through which 
the principal axes i and j are turned in the course of the shear is the angle a, so that 
after rotation i.e. in the final position, the angle between the rotated i-axes and the 
X-axis is 


B =2/4 — a/2 i) 
The components of strain with respect to the system x; or x, y, z can be derived from 


derived from Love’s equations, Art. 12, considering that ej = ej, = ex: vanish. 
Using the following equat.ons 


Cxxz = & cos2B + e; sin2B ) 
€yy =e sin2B + ej cos2B (7.6) 
Cxy = (e:—e7)/2 + sin 2B 


we find the expressions entered in the following Table VII. 


| 
| 
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TABLE VII 
Strain Components in Simple Shear 


H 


G 


—14V/24P22V 172 2VWVv+DiV—-L) 0 1-7/24+I2/2V/ 
2 ey = —1+-4//2—F2/2/ 0 -L2VinVtD(V—-L) —W/2 1-v/2—-I2/2vV 
Div T/2 WV inVvt+PMV-L) 2 TV 


We now go over to infinitesimal strain. Here one must be very careful. For instance 


what is —J°2 + W442 — 2, if I’ be infinitesimal? One cannot argue that with: 
infinitesimal I’ = y, y2 can be neglected against 4, so that the «/ expression becomes} 
2, and the whole reduced to —y?. The correct procedure yields from — J’? + 2(1+- 
+ [2/4)t—2 = —F2+ 2114+47°2/4 +...) —2 the expression —y2+ V 4 4 y2—. 
—2 —y2 + y2/4 = — 3y2/4. We thus find the expressions entered in Table VIII.. 

As can be seen, all measures result in the same first order expression for the shear- 
ing component of the infinitestinal strain tensor, equal to the one of the classical | 
theory. In addition, however, second order cross-terms appear. They are positive in} 
the x direction, the direction of displacement (vanishing in the case of the A-measure), 
and negative in the y-direction, the displacement-gradient (vanishing in the case of 
the G-measure). Note however that a positive strain component does not indicate} 
an extension, because in simple shear line-elements in the direction of the displace- 
ment (in our case x) do not change in length; nor does a negative strain compo- | 
nent indicate a compression for similar reassons. Such interpretations of the classical | 
theory are inapplicable. 


TABLE VIII 
Strain components in infinitesimal simple shear 


C G A A S 

1 be; y/2 

2-83 =z 

3 Ie y/4 y?/2 0 72/2 —y/4 
4 i, —¥2/4 

5 fea 3/gy? y2/2 y2/4 0 2/8 
6 Evy —1/gy2 0 —y2/4 —y2/2 3/gy2 
7 eey y/2 

8 & 0 


8. The very diverse expressions for strain components in different measures which | 
we have derived in the preceeding developments expressone single physical pheno- 
menon only, namely the strain ellipsoid with half axes 7;, 4;, Ay. Even the rotation | 
which is connected with simple shear is of no physical significance as long as we do. 
not define a coordinate system fixed in space by physical means. It is therefore diffi- _ 
cult to see what physical conditions are behind these expressions, especially when, | 
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as we have seen, there may be present a strain component, say ¢zx, without the line 
element dx changed in length through strain. This situation can be cleared up by 
proceeding to the loads and stresses causing the strains and connected with them in 
experiments by means of which one determines the parameters appearing in equa- 
tions of type (1.1) or (1.2). In these experiments specimens of certain simple shapes, 
mostly prisms or cylinders, are subjected to the action of certain simple loads. The 
mathematical problems connected with these experiments may be named primitive 
problems. 

The main primitive problems refer to homogeneous deformations, caused by (i) 
isotropic pressure or tension of a specimen of any shape, (ii) simple pull or push of a 
prism, (iii) simple shear. The classical solutions of these cases are 


(i) —p = XEy (8.1) 
which defines x, the bulk modulus, and 
(ii) o=Egys ae =—re (8.2) 


where o is the tensile traction, which defines E, Young’s modulus, and v, Poisson’s ratio, 
(iii) tT =y (8.3) 


where Tt is the shearing traction, which defines yw, the shear modulus. 
When E is found from the tensile test and ~ from simple shear x and y are calcu- 
lated from 


% = Eu/3(3u — E) 
ee ae i (8.4). 
We see that 
E = 2u + v) = 3x(1 — 27) (8.5). 


For mathematical convenience, but without physical signfiicance, Lamé’s constant 
expressed by 


%—2/3u=A, (8.6) 
is often used, as in Eq. (1.3). 


9. We now turn to the stress-strain relations. 
It has been shown by Reiner® that it the material is isotropic, i.e. the tensor function 
F of Eq. (1.1) is isotropic, it can be developed to 


Sim = Fo Oum = Fy Elm - Fy eja CAm (9.1) 


where the F are scalar functions of the invariants J, IJ; and //J,. Truesdell3 has 
expressed these functions as power series in the invariants, as e.g. 


Fa = Faooo + Fa1o0l + Fa200l2 + Faotoll + Fa30013 + Farol U + Fagoi Zl +... (9.2) 
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where a stands in turn for 0,1 and 2. 
Eq. (9.1) is thus replaced by 


sim = [( Foooo + Fo100!) 1m + Froooem | + [ ( Foz00l? + Fooiol!) drm + 
+ Fy,o9leum + Foooveia eam | + [ ( Fos00l3 + Foitol WI + Foot!) éim + 
+ (F 120012 + Fi0i011) Clim + F100! eja am | + eoeee (9.3) 


where the expressions within [ ] brackets are of the first, second and third order 
respectively, because the stress vanishes with the strain and the constant term F909 
is therefore identically equal to zero. 

We now go over to infinitesimal strain assuming that the F are all of the same 
order. 

If we compare Table VIII, we see that there are cases when J and ¢ are of second 
order infinitesimal, but J, JJ and JII will be at least of first, second and third order 
infinitesimal, ¢ at least of first order and ee at least of second order. 

In the case when Fojo0f dim + Fioo0€1m does not vanish in the first order terms, 
all other expressions in [ | brackets can be neglected and we have 


Stim = F100! bum + F000 €1m (9.4) 


as in classical theory. However, note that this may be the case for some definite ]and m, 
and then (9.4) is valid only for these /,m components. 
In the case when Fojo0/ 61m + Fio00 €1m vanishes in first order terms, we have 


Sim =(Fo100l + Fo200l? + FooiolZ) dim + (Fio00 + Fis00l) €1m + Fro00€ia Exam (9.5) 


However, when / itself is of second order, Foz99/2 can be neglected againstFo1 091 
and F;09l €im against F999 €im. We do not know of any problem when this is not 
the case, nor do we know of any problem in which terms of an order higher than 
two appear. In this connection it may be mentioned that the invariants J and J/ will 
both vanish only when there is no strain, or eim=0, in which case JJJ=0. We 
therefore omit the terms connected with the module Fo299 and Fj 100. 

We shall now deal with primitive experiments by means of wihch the moduli 
Fo100> Foo10, Fiooo and F999 can be determined. 


10. We first perform experiment (i) of §8, subjecting a specimen of the material to 
isotropic pressure -p, so that 


Sim = —P Om (10.1). 
Therefore from (9.4), to which (9.3) is reduced, considering Table IV, 


—P = 3y Foro0 + YF i000 = 3¥(Fo100 + F000 /3) (10.2) 
But from the same Table 


3y = Ey (10.3). 
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Therefore 
—P = (Fo100 + Fio00/3) €» (10.4). 
Comparing (8.1), we can therefore identify 
Fo100 + Fiooo/3 = K (10.5). 
11. The second experiment is (ii) the simple pull of §8, when 
00 0 
Sim =|| 0 0 0 (11.1) 
0 0cea 


Simple pull causes linear dilatation as defined in §6. We have from (9.3), which 
is reduced to (9.4), considering Table V, 


0 = Foio0 (Yt: —2yt) — Frooove \ (11.2) 
o = Foioo (v1 — 274) + Fiooovt 


Introducing 
gm=E; vit =v (11.3) 
we find 
Foio0 (1 — 2”) —Fiooov = 0 } (11.4) 
Foro0 (1 — 27) + Frooo = F 
From this 
Fio00 1 +») = £ (11.5), 
and comparing (8.5) , we can therefore identify 
F000 = 244 (11.6). 
Going back to (10.5) we now have 
Foi00 = K — 2u/3 = A, (11.7) 
and Equ. (9.4) becomes 
Sim = AI Sum + 2p em (11.8) 


as in classical elasticity, Eq. (1.3) 


12. We now proceed to simple shear as treated in §8. This is a case of plane 
strain when ézz, €yz and &zz vanish. The strain tensor is reduced to 


Exx Exy 


Elm = (12.1) 


Eyx Eyy 
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and therefore 


Egat + Exy” Exy(Exx +éyy) 
Ela Eam = (12.2) 
Exy(Exa + Eyy) Eyy? + Exy? 


From Table VIII this is reduced to 
y/2 Exe +Eyy 


EIA EAm = /2 (12:3); 


Exx + Eyy y/2 


We now calculate szy. From (9.5), (11.6) and Table VII we find 


Say = ylo:(2Qu +F 11002) + Faooo y/2° (€xa + Eyy) +... 12.4). | 


However, inspection of Table VIII shows us that expression ézz + €yy is of the 
second order in every measure. The second term on the right of (12.4) can therefore 
be neglected against the first term. Furthermore, Fj;99/ can be neglected against 2. 


The terms in the third bracket [ | of (9.3) can also be neglected, and so can all fol- | 


lowing | ]s. We therefore have in all measures 


Szy = PY (12.5) 


as in classical elasticity. 
It is different with the normal components. We find from (9.3), taking into con- 


sideration that J,ezzand eyy are of second order, and neglecting third and higher | 


order terms, 


Sra = ALI + 2m xx — Fooi0 y/4 + Froo0 y2/4 


Syy =A,I + 2m eyy — Fooio y2/4 + Fooo0 y/?4 (12.6) | 


Sez = A,I — Foo. y?/4 


In general, in order to produce simple shear, it will therefore be necessary to exert, 
in addition to the classical szy of (12.5), an isotropic “‘cross”-pressure of the amount 


Plo) = Fooi9 y2/4 — AT (12.7). | 


When the cross pressure p,,) is absent, the volume will expand, when acted upon 


by shearing tractions szy. This phenomenon was described by Reynolds in the | 
case of plastic deformations as dilatancy*. We may therefore term the factor Fooj0 | 


modulus of dilatancy (6). 
In addition ,there will be normal “cross’’-stresses sz.) and Syy(e) 


Syyle) = 2M Exe + Me y? \ 


Syy(c) = 2M Eyy + ue y? (12.8) 


* That such a phenomenon may be present in elasticity was suggested by Kelvin9. 
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where 
F000 = 4 be (12.9) 
and (Mc is the modulus of cross-elasticity. 


13. From the foregoing it appears that the general law of infinitesimal elasticity 
can be written as 


Stm = (AT + OT) dim + 2p Erm + 4Mtc E104 Ean (13.1) 


neglecting Fo299/? and Fyjo9/, and all third and higher order terms of Eq. 9.3). It 
is an open question whether there are cases where this approximation does not 
suffice. If they exist, they are not among those constituting the primitive problems. 
We have seen that in pure strain, second order terms are absent and classical theory 
reigns. However, when the strain is not pure, i.e. the principal axes are troated agianst 
a reference system formed by the external loads (x, y, z) second order terms are 
present even when the strain is infinitesimal. 

For the experimental determination of yu it is usual to employ an experiment, in 
which the external Joads are reduced to a torsional torque in the absence of any 
axial tension or pressure. As we now know, this will in general be accompanied 
by longitudinal and transversal positive or negative extensions. 

This can be verified through Poynting’s!° experiments. Poynting observed that 
cylinders of steel and copper in the form of long wires, when acted upon by forces 
equivalent to a torsional torque, are not only twisted as predicted by the classical 
theory of elasticity, but also elastically elongated and increased in volume. Reiner!! 
has treated Poynting’s observations in accordance with Reiner’s theory of finite 
elasticity using the Almansi-measure, and realizing only at the end of his investiga- 
tions that the strains in Poynting’s experiments were actualy infinitesimal. It is gen- 
erally believed that second order effects make their appearance in finite elasticity 
only, but we have shown in the preceding developments that this is not so. We shall 
treat Poynting’s observations in the light of our present investigation at another place. 


CONCLUSIONS 
It has been shown that when the strain of an elastic body is infinitesimal, i.e. when 
all components of the displacement-gradient tensor are infinitesimal, second order 
eflects will in general be present, in accordance with the stress-strain relation. 


Sim = (AI ar oll) Oim + 2u Elm + Atte Ela Exam 
with four elastic moduli A, 6, u and jc. The two moduli A and y are those of the clas- 


sical theory. In addition, however, two others are necessary for the complete de- 
scription, namely 6 the modulus of dilatancy and mw the modulus of cross-elasticity. 
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In infinitesimal pure strain the classical equation 
Sm = Al Oim + 2 Erm 


is sufficient. While 2 and uw can be determined by means of the customay experiments 
such as the tensile and torsion test the magnitude of 6 and uw: depends upon the mea- 
sure of strain used in defining the tensor e7m, of which €7m is the infinitesimal degenera- 


racy. 
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CONSTRAINED TORSION OF AN ELLIPTICAL CYLINDER BY FORCES 
APPLIED ON THE LATERAL SURFACE 


YEHUDA STAVSKY 
Division of Mechanics, Technion—Israel Institute of Technology, Haifa 


ABSTRACT 
St. Venant solved the free-torsion problem of a cylinder twisted by forces applied to the 
end sections. These classical boundary conditions occur seldom, if at all, in practical 
applications. A more practical case —constrained torsion— approximated by Foeppl. 
The free torsion of cylinders by tractions applied on the lateral surface was solved by Reiner, 
and explicit expressions for the stress in the case of an elliptical bar twisted by constant 
tangential tractions were derived by him. 

In the present paper the above solutions are superposed to include the case of constrained 
torsion of an elliptical bar fixed at one end (so that the cross-section remains plane at that 
end) and free at the other end from all tractions and hindrance, loaded by uniformly distri- 
buted tractions. 


1. St. Venant! (1855) was the first to solve the problem of “‘free”’ torsion of cylinders 
other than circular ones, where “‘free’’ refers to the unrestricted warping of all cross- 


sections. 
The displacements u; take the form 
Us; = — azy (1.1) 
Uy = azx (1.2) 
uz = ag(x,y) (1.3) 


where the torsion function g (x,y) is harmonic throughout the cross-section R of 
the twisted cylinder, while its normal derivative assumes a given value on the boun- 
dary; a is the twist per unit length and z is the axis of the cylinder. 

St. Venant’s solution requires that: 

(a) the surface of the cylinder be entirely free from loads; 

(b) the tangential tractions szz, Szy be distributed in a specific manner, which 
is the same in any cross-section ; 

(c) the only non-vanishing components of the loads be reduced to the torque Mz. 

In practical applications the assumption that the torsion is free, seldom, if ever, 
corresponds to the actual physical situation when, at least, one end of the cylinder 
is fixed and the torsion is thus “constrained”. 

An extension due to A. Foeppl? (1920) imposes on St. Venant’s solution the re- 
quirement that one of the sections of the twisted cylinder remain plane, L.€..Uz| ge0== 0. 
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Following the principle of St. Venant — “le mode d’application et de réparti- 
tion des forces vers les extrémités des prismes est indifférent” — Foeppl assumed that 
there is a local effect due to the constraint in the section z = 0, so that the axial 


displacement is 
uz = a.9(x,y)(1—e7*”) (1.4) 


where f is a constant. Using energy methods, he approximated this problem for a 
beam of elliptical cross-section (Figure 1). His expressions comply with the equa- 
tions of equilibrium. They do not, however, all comply with the equations of compa- 
tibility. 


——- wee > Xx 


Figure 1 


In all the foregoing considerations the lateral surfaces of the cylinders are free from 
external load. 

M. Reiner 3:45 (1915, 1925, 1933) first investigated the free torsion problem of a 
cylinder twisted by forces applied on the lateral surface. 

The given cylinder has an arbitrary cross-section R and the boundary curve C, 
the unit normal to C being denoted by ». 

Let the system of coordinates be chosen as shown in Figure 2, with the fixed 
end of the cylinder in the plane z = 0, and the free end in the plane z = 1. 


y y 
\ 1, 


4 


Figure 2 
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The prescribed surface tractions S; (i = x,y), assumed as independent of z, are 
_ directed parallel to the xy-plane, and produce a twisting moment m, per unit length. 
Thus, the boundary conditions are: 


Sz=0, Jf,S(x,y)ds =0, f {(xSy—ySz)ds = mz (1.5) 


The relaxed boundary conditions at the ends take the form: 
(a) at the end z=0 


SJsacdA = [fsyxdA = ff szedA =0 


R R 


. 
SS@xsyz — YSzx)dA =m;z:1= M; 
; (1.6) 


Sf xs22dA == SfvszedA == 0 


R R 


(b) at the end z =/ 


SJ seedA = [J xse2dA = ff yszedA =0 
R R R (1.7) 


Szz ==1((). Szy =="() 


The boundary conditions on the lateral surface are 


(1.8) 


Sx COS(X,”)+syz cos(y,v) = Sx(x,y) | 
Sy COS(X,) +Syy Cos(y,v) = S(x,y) 
Sze COS(X,v)+sSzy cos(y,v) = 0 
The problem is to determine a system of tractions s;; which satisfies the equilibrium 
equations (the body forces F; are assumed to vanish) 


Sigg =0, (Gj = 1,2,3) (1.9) 


the relaxed boundary conditions (1.6), (1.7), the boundary conditions (1.8), and the 
Beltrami-Michell compatibility equations. 


1 
7214 es 6, 4 =0 (1.10) 


where Sj; is the stress-tensor, 6 its first invariant and o Poisson’s ratio. 

Reiner uses the St. Venant semi-inverse method of solution, which consists in 
making certain assumptions about the stress distribution, leaving sufficient freedom 
to permit the satisfaction of the differential equations and the boundary conditions. 

It is obvious that the torsional couple transmitted across each cross-section is a 
linear function of the distance from that cross-section to the free end of the cylinder; 
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hence it is assumed that the stresses in the cylinder have, in part, an pet | 
similar to those in St. Venant’s torsion problem. 

As the applied tractions do not vary along the length of the cylinder, it is dentro 
that the deformation is independent of z, or 


Uz =Ux(x,y), Uy =Uy(x,y), uz =0 


For such deformation, known as plane strain, it is necessary that szz = 0, but as an 
variable sz, at the end sections may produce a bending or stretching of the cylinder! 
we have to add the bending and torsion type solutions. 

The 543 system obtained is 


Sa Po \ 
a’ ae D Oy2 
ee a ay 
a BP) ax2 
a = 0 py —p + A+Bxt Cy | 
st x (1.11)| 
fone. eee 
ae = HI (Se + 3) 
fe ig 
Set = 4 —2(5% — ») 
a’ d(x y ) + oxdoy 
where y7y(%y) =0 in R 
dp 
7 = 7 c0s(x,») —x cosy”) on C (1. 2 | 
yyp=0 inR | 
d a yeaa ty ) 
ae = ae cos(x,”) (a T >) cos(y,r)— 
ds ‘ dy 2 4 a c | 
d , &y x2 —y2 S Y « ee | 
Tay oe), aa ei ee ee 
S (xSy — ySz) ds mz 


Qaa=--__eeo—-: RrwmRRTTT- 2 
5 P ; (1.14) 
4JI (Pty tae _ ys )dxdy 


Thus, if the torsion function p is known for R, the value of the derivatives of y is 
known on C, and the problem is reduced to the fundamental boundary-value pro- 
blem in the biharmonic equation. 
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2. Explicit formulae are known only for the stress in a beam of elliptical cross- 
section twisted by constant tangential load S. 

The bi-harmonic function y (x,y) is 
b2xy3 —a2x3y 


y(x,y) = 12(@2 + By (24) 


where a and b are the semi-axes of the ellipse: The torsion function @ (x, y) is 
b2—q2 


(x,y) = mip ? (2.2) 


The stress components, following Reiner, take the form 


Seer be ee 7 
Szz = 2(a2-+62) y 
a a’ h2 e 
Syy = 2(a2+-b2) y 
_ — F(a? (o+2)_ 
Ses = a+b)” (2.3) 
'h2 
Syz = re —z) 
a’ a2 
Szz = = pap 2) 
2) ee ee ee 
Say 4(a2-+b2) (b x a y?) 
a2h2 7 
where a ‘= Ka2b) Mei 2S52ab 


An extension of this theory to the problem of constrained torsion is given here for 
the case of an elliptical section twisted by constant tangential tractions s applied on 
the lateral surface. Then, making use of the linear relation 


&ig = 4(ui,g + 5,2) (2.4) 


we find by integration the following expressions for the displacement components 


te = ESTP (4a +0) +08? ]x4y-+(@2 +2) 9(1—2)°— 
— [244 [62(1 +0) —oa2J]y3} 
uy = ee { § [b2(1 +0) —oa? Jxy2—(a2-++b2)x(I—z)2 +r 2-9) 


+ $[02—4 [—a%(1 +0) + ob? }]x3} 


‘1+o) 62—a? 
a Sig en “grnpooe *) J 


132 YEHUDA STAVSKY Bull. Res. Counc. of Israel! 


where the constants of integration representing a rigid body translation and a rigid 
body rotation can be assumed to vanish by fixing the origin, an element of the z-axis 


and an element of the xz-plane at the origin, These give six independent conditions 
Oux Uy ou 


x Pie een 6) 
ee 0 at (0,0, 0) (2.6) 


Ug = Uy = Uz = 


From Equ. (2.5,3) it can be seen that the longitudinal displacement uz varies linearly 
along the cylinder and becomes a maximum at the fixed end (x = 0) in the same 
manner as the torque (Compare Figure 3). 


Figure 3 


U, diagram for xy = const. 

3. The problem of constrained torsion by forces applied on the lateral surface: 
can now be solved with the help of the three following problems: 

(a) the free torsion, where the bases of the cylinders are subjected to external] 
tractions, following St. Venant; 

(b) the free torsion of the cylinder by forces on the lateral surface, following Reiner; / 

. (c) the constrained torsion with the lateral surface free from load. 

While the solution of (a) is known for many sections, the solution of (b) is only) 
known for an elliptical beam twisted by constant tangential tractions. The problem) 
(c) was only approximated; by Foeppl for the elliptical beam, by S. Timoshenko’! 
(1922) for a narrow rectangle and by J. Nowinski8 (1954) for a rectangular bar. | 

Let us construct the solution for the constrained torson of an elliptical bar twisted! 
by a constant tangential load applied on the lateral surface. 

We impose on (b) a moment M; according to (a), in order to satisfy the boundary) 
conditions uz = 0 at z = 0. Its magnitude is 


4S(a2+62 1 
Dm Dot Fh shay) 1=2nablS (3.1) 


Me = D "0.— Da, E bE E 
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where D, the torsional rigidity of the elliptical cross-section, has the value 


mua3b3 
and uw is the modulus of rigidity. 
Using Foeppl’s approximate solution we comply with the relaxed boundary condi- 
tion at z =/ (1.7), not interfering with the already satisfied boundary condition at 
-z=0 (uz = 0), and so the problem (c) is solved. 


The szz system takes the form 


Fe PUES Lee 
ie a’ b2 

Syy = Q(a2-+-b2) xy (3.4) 
a’ (2+0)(b2—a?2) 

Sez. = i ee xy + a y(1+o)! ——— aa aye" C= (aes) 
a’ b2 1 

Syz = a2+b2 x(J—z) ae e~Pa, 4Kb2x(a2b2 —52x2= a2y2)0 (3.6) 
a’ a2 8 2 

Sez = a ee Pa ae 2. 4Ka2y(a2b?—b2x2—ay2) (3.7) 


Saye (ELS (b2x2—q2y2) + K(a2b2 —h2x2 —q2y2)2¢e— Bz (3.8) 


where 
) b2—a2 1 
K = a@f3(1+0) 1_— 16(a Eb) as (3.9) 
| 2452 (a2+62)2 16 
2 = a a. ——— 
abp te ee (3.10) 
: 4S(a2+b?) 


4. Equs. (3.3) to (3.11) give expressions for the stress-components in a cylinder of 
‘elliptical cross-section, fixed at the end z = 0 (so that the cross-section remains plane) 
and free from all tractions and hindrance at the end z =/, loaded by tangential 

tractions S applied on the sides of the cylinder in planes parallel to the xy-plane 
and constant along the circumference of the ellipse and along z. The resultant of 
this loading is equivalent to a total torque Mz = 2zabls, which is taken up by trac- 
tions acting tangentially in the cross-section z = 0, and distributed in accordance with 
(3.6) and (3.7). 

When we introduce a = b = R we find Reiner’s solution for a circular bar where 

warping does not occur and no constraint is therefore required to keep the section 


plane. 
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Five problems may be indicated for further study: 
(a) explicit formulas for Reiner’s problem in the case of sections other than a 
ellipse; 
(b) exact solutions for the constrained torsion problem approximated by Foeppl! 
Timoshenko and Nowinski; 
(c) constrained torsion of cylindrical bars by constant forces on the lateral sur- 
face in the case of an arbitrary cross-section; 
(d) extending (a) to the case where the applied surface tractions vary along the 
length of the cylinder; ' 
(e) extending (c) in accordance with (d). 
This paper was submitted to the Technion—Israel Institute of Technology ag 
a thesis, in partial fulfilment of the requirements for the Master of Science degree. 
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SAFETY FACTOR IN STRUCTURES IN THE ELASTIC RANGE 
WITHOUT STRESS-LOAD PROPORTIONALITY 
A. ZASLAVSKY 
Division of Structural Engineering, 
Technion-Israel Institute of Technology, Haifa 


| ABSTRACT 


The paper deals with elastic design of structures having no stress-load proportionality. It 
I is suggested that these structures be designed on the basis of the load producing a critical 
} stress, rather than on the conventional basis of the admissible stresses. The suggested 
method ensures the intended required safety factor against the appearance of the critical 
stress (for instance, yield stress). 


GENERAL 
| This paper deals with the design of a group of structures by the conventional 
j elastic method, i.e. on the basis of admissible (working) stresses. Ultimate load de- 
| sign is mentioned only for the purpose of comparison. 

| In these structures there exists no proportionality between load and stress, in spite 
| of the stress being within the elastic range and the validity of Hooke’s law. Examples 
_ of such structures are: structures where the change in their geometry due to deform- 
| ation cannot be neglected in the application of the equilibrium conditions; composite 
| structures; partially precast structures ; strengthened structures (beams, columns, found- 
| ations); structures erected by stages; partially and totally prestressed structures; struc- 
| tures with thermal and assembly stresses; structures with “conditional restraint”, etc. 
It is maintained that the design of such structures should be based on the failure 
) load producing a critical stress, rather than on the admissible stress; only this method 


| ensures the intended required safety factor against the appearance of the critical 
| stress (for instance, yield stress) which is the true strength criterion of the elastic 
| design. It should also be noted that in the above structures only a modified superpo- 
) sition method can be applied. 

_ According to the conventional elastic method, the admissible load Paam is the 
- one producing the admissible stress oaam in the critical section of the structure. 
' Through long routine use, the concept of admissible stress has come to be re- 
, garded as an independent criterion in its own right; in fact, however, the admissible 
| stress depends upon two other quantities and equals the critical stress ocr divided 


| by the safety factor k: 


Oadm = ar/k (1) 


| Received August 8, 1957 
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ohn 
We shall not deal here with the numerous problems related to these two quanti] 


ties. There is no agreed definition of the critical stress. In steel structures the term) 
usually refers to the yield stress oy; in concrete, it is mostly the prism (or cylinder} 
strength*. The safety factor is governed by a number of factors, but will be usec 
here as defined in existing elastic Codes of Practice — or by definition (1), ocr ancy 
Oaam being given by the Codes. | 

For simplicity, assume that the structural steel has a yield stress (= proportional 
ity limit) of cy = 2400 kg/cm? and oaam = 1200 kg/cm?, k = 2 being the safety factor? 
By the conventional method, Paam would then be the load producing oaam = 
+ 1200 kg/cm? in the steel structure. On the other hand, the failure load P, is the one 
producing the yield stress oy = + 2400 kg/cm? at the same point; hence the 
admissible load Paam (the asterisk added for distinction from the conventiona} 
Paam) should be: | 


Paam = Py/k (23) 


In the ordinary cases, with o proportional to P, Paam = Paam; but when no sucht 
proportionality exists, Paam S S Paam. In these cases it is suggested that the design 
be based on Pagm. This method is illustrated by a number of examples. i} 

This paper does not deal with cases where the equivalent stress is determined im 
accordance with the different strength theories. It is also assumed that problems of} 
stress concentration and fatigue do not arise, and that limited dynamic loading is 


expressed by an equivalent statical loading. 


EXAMPLES 


Figures 1 and 2 show two well-known cases, where equilibrium conditions should 
be applied to the deformed structure. Figure 1a shows two horizontal bars connect- 
ed by hinges and supporting a vertical load P. The load-deflection curve is concave, 
and so is the load-stress curve, given in Figure 1b (for A = 1 cm2 =const; E=2.1 x‘ 
10° kg/em2): | 


tI can be seen that: 


The existing (working) stress produced by Paam is o* = 1510 kg/cm2 > ogam, butt 
a safety factor of 2 against the appearance of oy is still ensured; for such a reliable: 
material as steel, even higher working stresses should raise no objection. 

Figure 2a shows a slender bar subjected to two equal and opposite eccentric forces} 
P. Here the secant formula applies (Figure 2b): 


* In buckling, it is the critical buckling stress. 


Vol. 6C, 1958 A. ZASLAVSKY 137 


Rage oo 5s ae ne: 


y ie 

AS (a) 

| B=230 
+ 
(eta =475 
Baim = 4 (b) 
(ORG kgk im) 
| S- eNS 
| NY + 
| > N 
| ‘s u 
| Cae © 
. Figure 1 
() 
| 
oe 
Cache, oy spies: 

| Figure 2 
| p°\ Peséc [4/P/(ED] 
| Omaz = =~ + fe (4) 
A Z 
) The curve is convex and design by oaam (or Paam) would not provide an adequate 
| safety factor against oy. 
| A simple model of a prestressed structure (Figure 3a) consists of a short pipe 
| (A, = 1 cm?) with a bolt (A, = A; = 1 cm?) passing through it. By turning the nut, 
| 
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(a) 


(5) 


Figure 3 


the bolt is pretensioned (o,, = 1000 kg/cm?) and the pipe precompressed (o,,, = | 
—o;,, = — 1000 kg/cm?). We now determine the admissible axial tensile force P — 
which may be applied to the pipe. Figure 3b shows the P/o, and P/a,, lines* : | 


Paam = 1400 kg > Paam = 400 kg 
We could also have made our point, of course, by prestressing up to o,, = — o,,, = 
1200 kg/cm2 = ogam; in that case: Paam = 0 but Paam = 1200 kg. 


The dashed lines in Figure 3b show the extension into the plastic range. The 
ultimate load Per of the structure does not depend on the prestress: 


Por = Ayoy + Ayoy =2X 1.0 x 2400 = 4800 kg 


* Up to the point where the pipe’s prestress is reduced to zero (P=200 kg), the force P could have 
been applied directly to the bolt without affecting the P/o lines. 
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A composite steel-concrete beam is shown in Figure 4a. This is one of the longi- 
tudinal beams supporting a bridge with a concrete roadway. The idea is to use the 
steel beam (I = 813,000 cm4) for carrying the dead load in the first stage of erection 
(bending moment M,). In the second stage, the composite beam (‘transformed 
section”; Ieq = 2,046,000 cm4) carries the additional live load including pavement 
etc. (Mp). Special welded dowels ensure the connection between the concrete slab 
and the steel section. The lower flange of the steel section is made larger 3 in order to 
achieve better utilisation of both flanges. 

Assuming the bending moment of the first stage to be My = 100 mt, we determine 
the additional admissible bending moment My which may be carried by the compo- 
site beam at the second stage. 

The stress distributions which lead to Mpaam) and to Mgaam) = Mg+ Mp aamy 
by the conventional method are shown in Figure 4a*. 

The M/o lines for the extreme fibres of the steel section are shown in Figure 4b, 
which also reveals: 


Ma adam) = My/k = 387/2 = 193.5 mt > Maam = 152mt 
Magus = M qaam) — M,= 93.5 mt > Mpaam) = 52 mt 


The (working) stresses produced by Maile ) = 193.5 mt are o, = 1310 kg/cm? 
and 07 = 1280 kg/cm?. 

When actual loads are expected to exceed design loads, strengthening of existing 
structures may become necessary. Examples are: steel joist with additional welded 
flange-plates; truss-strengthened beam; arch with tie-rod, etc. 

Figure 5 shows a symmetrical model with “conditional restraint” consisting of 
three steel bars (A = 1 cm? = const.) supporting a vertica] load P. The bar X comes 
into action only after contact has been established at B (P = P;). The additional 
load P—P; is then resisted by all three bars, X resisting twice as much as Y (for the 
given slope of 45°). 

Figure 5 also shows the relationships P/o, and P/oy. The set-off continuous lines 
represent the case t =0, the set-off dashed lines showing the extension into the 
plastic range. The change of ¢ permits the regulation of the stresses in the bars. In 
this example it is possible to utilise simultaneously the admissible stress in both 
bars when ¢ corresponds to an initial stress of oy, = 600 kg/cm? (continuous 
lines); in this case we arrive at: 


Paam = 4940/2 = 2470 kg < Paam = 2895 kg 


If, on the other hand, ¢ is made to produce oy, = 800 kg/cm? (point b; dot-dash 
lines), Paam is governed by bar X, while Paam is determined by bar Y: 


Paam = 5230/2 = 2615 kg > Paam = 2495 kg 


* Since the stresses in the concrete slab are quite low, we are justified in neglecting the creep of the 
concrete (for M,) and also the non-proportional stress-strain relationship. 
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A structure with conditional restraint can also be considered as partially prestressed 
(bar Y having initial stress before X comes into action). Another example: a pipe (sub- 
jected to internal radial pressure) inside a non-tightly fitting sleeve. Partial prestressing 
is especially used in bridges. Consider, for instance, a three-span beam (Figure 6), 
the middle span BC being usually much longer that the other two. In order to re- 
duce the moment at m it is advantageous to let the end spans AB and CD act as 
freely overhanging cantilevers in the first stage and establish contact with the outer 


supports A and D in the second stage only. 


A LA 
— 


Figure 6 
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The simpler case of a cantilever with a conditional support B is illustrated in 
Figure 7. The P/M lines also represent the P/omaz lines up to M = My. The dashed 
lines show the extension into the plastic range, neglecting the curvature of the 
elasto-plastic stage. 

CONCLUSION 

Non-proportionality between stress and load in the elastic range has been 
illustrated in a number of examples. 

The relation may be represented by graphs of ditferent shapes — curve (concave 
or convex), straight line not passing through origin, broken straight line, etc. In the 
majority of cases design by Paam == Py/k entails considerable economy, but there 
are also opposite cases. Therefore each case should be studied in order to ensure 
the required safety factor. 


EXPERIMENTS ON AERATION METHODS 
DURING MYCOLOGICAL CITRIC ACID FERMENTATION 


Z. BERK AND J. B. S. BRAVERMAN 
Division of Food and Biotechnology, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 
In order to create better contact between air and the fermenting media during the produc- 
tion of citric acid by submerged fermentation, three different methods of aeration were tried: 
a) Fermentation on artificially-produced foam. 
b) Fermenting while spraying the medium into air. 
c) Aeration of the fermenting medium with the aid of a special vibro-mixer. 
Special fermentors were assembled in each case and series of runs were tried with sucrose 
solutions and diluted molasses. 
Of the three methods only the last one, using the vibro-mixer, gave positive results, 
reducing the time of fermentation and achieving quite a high rate of conversion. 


INTRODUCTION 

The first trials in producing citric acid by the mycological method, although 
unsuccessful at the time, actually began 65 years ago following the publication of 
Carl Wehmer’s! first observations, namely, that certain microorganisms isolated 
from the air produced this acid from sugar. But only in the late thirties of the present 
century this method became an important factor in the manufacture of citric acid. 
To-day it has practically replaced the old Scheele method of producing citric acid 
from lemon juice, and in the U. S. A. alone the annual output of the biological 
product has reached a total of over 25,000 tons. 

Important advances have accompanied the development of this industry. These 
have been ably summarised in a number of articles23. The original species of 
microorganisms described by Wehmer were replaced by more productive fungi, 
such as Aspergillus niger4 and Asp. wentii>. Fermentation was originally, and still 
is in many cases, carried out in shallow pans® and only recently various experiments 
in submerged fermentation were attempetd?.8. As the development of these fungi 
requires an abundant supply of air, and as ordinary aeration procedures were inade- 
quate, the shallow-pan method was usually preferred in spite of the fact that it 
entails a large initial outlay and excessive maintenance charges. This method of 
producing citric acid is also comparatively slow (8 to 10 days) and is, therefore, 
open to danger of contamination by undesirable microorganisms. 

Many attemps have been made to provide the necessary aeration by such means 
as increasing the surface exposed to the air by soaking wood shavings, or similar 
porous material, in the medium to be fermented? 


i Febru 9, 1958. 
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In the present work three attemps have been made to overcome the difficulties 
of aeration, and although they were not entirely successful it was considered worth- 
while to describe them herewith. Studies of submerged fermentation of citric acid 
published to date are all based on one of two methods, namely the shaking flask 
method or the column-fermentor, or on modifications of both. The object of the 
present work, however, was to study three other possibilities: 

a) Fermentation on artificially produced foam. 

b) Fermenting while spraying the medium into the air. 

c) Aeration of the fermenting medium with the aid of a special vibro-mixer. 

The working theory underlying all three of these methods was to permit better 
contact between air and the fermenting media, thus reducing the period of active © 
fermentation and preventing the development of undesirable microorganisms. 


EXPERIMENTAL 


A. Fermentation on foam 


Cahn9 reported considerably accelerated rates of fermentation by Aspergillus 
niger grown on solid materials, such as wood shavings or beet pulp, impregnated © 
with the medium. Such accelerated rates were no doubt due to greater contact surfa- | 
ces between the three phases: fungi, medium and air. The final products of fermenta- 
tion were then extracted by diffusion. On the other hand, it is well known that liquids | 
used as media in this process are apt to produce considerable quantities of foam 
during aeration. In industry such foam is combated by various means, such as 
addition of anti-foam or mechanical foam-breakers. 

The idea was to utilise the foam as a bed for this fermentation, and grow the 
fungi on the enormous surface created by it, thus dispensing with the addition of a 
solid phase and consequently also eliminating the necessity of extracting the citric 
acid. 

The apparatus used (Figure 1) consisted of a large fermenting vessel A, only partial- 
ly filled with fermenting medium. The air entering through B facilitated the formation 
of copious foam which was normally trapped in vessel D, broken up mechanically 
and returned to the main vessel while the air escaped through F. 

In all experiments by this method the medium consisted of: 160 g sucrose, 1.0 g 
KH POq, 2.5g NH4NO3,0.25g MgSO,4. 7H20, the whole made up to 1 litre with distilled 
water. The medium was sterilised at 15 psig for 15 minutes and the acidity adjusted 
to pH 4.0 with HCl. As foaming agent, white of egg was used. Before each 
run the whole apparatus was sterilised with a 10% solution of formalin, washed 
with hot water and detergent, rinsed several times with hot water and finally twice 
with sterile distilled water. 

Inoculum used: 20ml of a heavy suspension of spores of the strain Aspergillus 
niger ATCC 1015 cultivated on a solid medium of sugar-salt-agar similar to Czapek’s 
medium, four transfers having been made before final sporulation. 
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Figure 1. 
Foam fermentor assembly 


A - foam chamber 
B - air inlet 

C - foam run-over 
D - foam breaker 

E - motor 

F - air outlet. 


After inoculation, aeration was started at a rate of 100 ml per min. Almost imme- 
diately a permanent foam built up and continued to circulate slowly. Samples were 
drawn every 24 hours and tested microscopically as well as for total acidity and 
residual sugar. 

Four identical runs, six days each, gave negative results as shown in Table I, 
which represents one typical run. Practically no acid was formed and even the growth 
of fungi was very poor: formless flakes were found in the liquid below the foam. 


TABLE I 
Time Acidity Appearance Residual sugar 
hrs. % citric acid under microscope % sucrose 
0 _ clear 14.0 
24 oa 2 14.0 
48 sess ss 14.0 
72 _ Germination sets in 14.0 
96 0.3 Non-uniform mycelial 13.0 
flakes 
120 0.3 oS a 13.0 


144 0.6 i 12.0 
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The failure of the above experiments was explained by the possibility that the 
presence of the foam impedes the development of the mycelium; it was assumed 
that with the use of previously prepared mycelium. added to the medium in finely 
disintegrated form, better results might be obtained. 


The mycelium was then prepared as follows: 

100 ml of the sugar-salt medium were inoculated with a suspension of spores prepared 
as described above. After 5 days’ incubation at 28°C, a continuous compact mycelium 
was formed without sporulation. The liquid was decanted and the mycelium, washed 
with 3 portions of sterile distilled water, was transferred to a Waring Blendor together 
with the substrate intended for the fermentation. Disintegration produced a fine 


suspension of the mycelium in the medium. This suspension was used as inoculum > 


to 1 litre of substrate, all conditions of this experiment being the same as above. 
As control, 100 ml of the substrate were inoculated with a small portion of the 
disintegrated mycelium and left for incubation at rest. However, while in the stationa- 


ry control flask a new continuous mycelium was formed and after 6 days 4.4 g of © 
citric acid were present, the main experimental batch again proved unsuccessful: | 
the disintegrated mycelium particles agglomerated into formless flakes and failed | 


to rise with the foam. Practically no sugar was used up by the fungi and no acid 
formed. A duplicate test under similar conditions was equally unsuccessful. 


B. Fermenting while spraying the medium into air 


This experiment was based on the assumption that better results might be obtained 
by spraying the medium into air, as is done, for instance, during hydrogenation 
of oils, in which the liquid phase is sprayed into the gaseous phase. 

Figure 2 shows the specially-designed sprayer, as used in the apparatus shown 
in detail in Figure 3. The liquid in flask A was drawn by sprayer B, through which 
a large volume of air was passed, distributing the medium in the form of a fine 
mist of droplets. The mist-carrying air escaped through a series of cyclone-shaped 
traps (E) in which the droplets were separated and returned to the fermenting vessel. 
At the same time fresh air was introduced through pipe C in order to agitate the 
liquid. 

The following procedure was adopted: 

One litre of the sterilised sugar-salt medium was transferred to vessel A. After inocu- 
lation with a suspension of spores, aeration through pipe C began at the rate of 
500 ml per minute. Aeration was maintained for 48 hours, when germination set 
in and small particles of mycelium were formed, giving the entire liquid an appearance 
of a cloudy suspension. At this point spraying was started, with a discharge rate 
of at least 21 litres of air per minute. Sterile water was constantly added, to make 
up losses due to evaporation. 


The results obtained in 5 runs are given in Table II. 
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Figure 2. Figure 3. 
Sprayer Spray-fermentor assembly 
A - spray chamber B - sprayer C - agitating air 
D - air for spraying E - traps F - air outlet 


C. Aeration through vibro-mixing 


The use of a vibro-mixer, practised in chemical laboratories, was recently adopted 
also in biotechnological work. Such vibro-mixers are manufactured at present in 
several industrial sizes. In these experiments, a laboratory 40 W model of Swiss 
make was used. This mixer contains no rotating parts; instead, a large iron armature 
is attracted by an electromagnet supported by springs (see Figure 4), giving rise 
to rapid vibrations in synchronism with the electric current (in this case, 50 cycles 
per second), the amplitude of the vibrations being too small to be observed visually. 
When the vibro-mixer is used in conjunction with a special aeration tube, the entering 
air is finely dispersed in the form of minute bubbles. As far as known, such vibrators 
have not yet been used in citric acid fermentation. Figure 5 shows the apparatus 
used in the following runs. 

Before starting the experiments on molasses, trial runs were made on model 
sugar-salt solutions containing, per litre, 170 g sugar, 2.5 g NH4NO3 and sufficient 
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Time Acidity Residual sugar Aeration Remarks 
hrs. % citric acid % sucrose ml/min 
Run No. 1 
24 — — 500 Germination sets in 
48 — — 500-+21,000 Small units of mycelium form 
72 0.22 12.7 500+21,000 
96 0.64 11.2 500+21,000 Mycelium agglomerates into 
large clumps 
120 0.70 11.2 500+21,000 System clogged 
Run No. 2 
24 — _ 500 Germination sets in 
48 — oo 1000+21,000 Mycelium develops 
72 0.30 12.0 1000+21,000 Mycelium develops 
96 0.98 11.7 1000+21,000 Mycelium develops 
120 1.70 10.9 1000+21,000 Mycelium develops 
144 2.80 9.3 1000+21,000 System clogged 
Run No. 3 
24 — — 500 Germination sets in 
48 — — 1000+21,000 Mycelium develops 
72 0.18 13.3 1000+21,000 Mycelium develops 
96 0.20 13.0 1000+21,000 Mycelium develops 
120 0.18 12.4 1000+21,000 System clogged 
Run No. 4 
24 — — 500 No mycelium formed, 
48 — = 500 liquid remains clear 
144 — —_ 500 
Run No. 5 
24 oo — 1000 Germination sets in 
48 — oe 1000+21,000 Mycelium develops 
72 0.32 12.0 1000+-21,000 
96 0.83 11.6 1000+21,000 
120 1.05 11.0 1000+ 21,000 
144 2.8 9.3 1000+ 21,000 
168 2.8 7.0 1000+21,000 System clogged, strong odour 
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of alcohol 


HCI to bring the solution to pH 2.0. Aeration and mixing were started immediately 
upon transfer of the hot solution into the fermentor. After cooling, the solution 
was inoculated as usual and the aeration rate brought up to 500 ml per min. Acidity 
was measured once a day. Fermentation ended after 5 days, with a yield equivalent 
to 48% conversion. 


In the last stage of these experiments, sugar-beet molasses were used as a source 
of carbohydrates. Their composition was as follows (Table III): 
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Figure 4. 


Vibro-mixer 
(A. G. fiir Chemie-Apparatebau, Ziirich) 
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| TABLE II 

| Specific gravity 1.455 “H 
| Water 17.5% 
Total solubles 78.7% 

_ Non-solubles 3.8%, 

; 

i ‘otal sugars SOYA 
Ash 6.1% 
 Fe,03 0.0488 % 
20s 0.0894 % 

| Figure 5. rd 


_ Fermentor assembly with vibro-mixer 


A - fermenting vessel, B - vibro-mixer 
'C-air inlet, D - diaphragm seal 
| E - anti-foam receptacle F - sampler 
G-constant-temperature bath 

(all necks cotton-plugged). 
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Treatment of the molasses before fermentation consisted in their dilution to 
22° Brix (approximately 15°% sugars), sterilisation at 120°C for 20 min. and preci-- 
pitation with potassium ferrocyanide (0.6 g/1). 

The diluted molasses were used for propagation as well as for the citric acid fermen-- 
tation. For propagation purposes, 100 ml batches of diluted molasses were seeded! 
with 1 ml of a suspension of 4-5 days old spores after fourth transfer to sugar-salt- 
agar slants. Growth was achieved in flasks, by the shaking method, at 28°C. After: 
2-3 days very small pellets (0.05-0.1 mm) were obtained. 

The entire mycelium, together with the solution in which it had been obtained,, 
was used as inoculum for | litre batches of diluted molasses in the fermentor. Before 
inoculation, the solution was aerated at increasing rates, from 10 to 500 ml per 
min. per litre. This stage was accompanied by considerable foaming, and the rate o 
aeration was increased further after the foam had subsided each time. Although: 
this preliminary aeration helped to check the formation of foam during the subse-: 
quent stage, it was still found necessary to add, after inoculation, about 10 ml of 
a 2% solution of stearyl alcoholin paraffin oil to prevent excessive foaming durin; 
the first ten hours of fermentation. During the fermentation stage, the rate of aeration 
was 500 to 600 ml per min., at a constant temperature of 28°C. 

Samples for analysis were drawn every 24 hours, with the following results. 


TABLE IV 
Run Dura- Initial Final Acidity Weight Con- Yield, based | 
No. tion of sugar sugar VA of ver- on total | 
fermen- concen- concen- citricacid mycel- sion sugar 
tation tration tration ium 
days Yo % g/100 ml Ye 
1 5 15 33 9.80 2.0 83 
2 6 15 3:3 10.00 1.3 91 
3 5 15 3.8 10.5 1.5 93:5 


DISCUSSION OF RESULTS 


1. The foam method | 


Contrary to expectations, Aspergillus niger did not develop on foam. This can prod 
bably be explained now as due to one or both of the following causes: 
a) This fungus is not a monocellular microorganism, and in order to be capabl 
of metabolic activity the mycelium should be composed of a certain minimum num- 
ber of cells. In other words, there should probably be a certain minimum thicknes 
of the mycelium unit (partly immersed in liquid) which the foam film is too thi 

to support. 
b) Surface tension forces may also play an important role here, and are highl 
likely to preclude any metabolic exchange between the cell and the foam film. 
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2. The spray method 


_ This was actually abandoned for technical reasons without conclusive results. In 


the first place, it was found impossible to prevent the formation of agglomerations 
of mycelium (due to insufficient agitation of the liquid by air alone) which clogged 
the system after four to six days. 

Secondly, the spraying process itself required such large amounts of air, that some 
of the runs failed owing to excessive pressure in the system. 

The results of runs 1, 2 and 5 show, however, that the method is feasible in prin- 


_ ciple, although the system is somewhat cumbersome and requires a fermentor with 


j 


a capacity of 15 to 20 times the volume of the substrate. Since fermentation by this 


_ method could not be accelerated as expected, it was finally abandoned in favour 


of the vibro-mixer. 


3. The vibro-mixer method 


The vibro-mixer proved extremely effective and gave quite satisfactory results. Both 
when using sugar solutions as well as molasses, the results compared favourably 
with those reported in literature, and experiments on a larger scale should, in the 
opinion of the authors, prove rewarding. The time of fermentation was reduced 


to 5 days as against 10-11 days by the shallow pan method, and the rate of conver- 


sion was high. The only difficulty encountered was in the formation of foam immedia- 
tely after inoculation, but this can no doubt be avoided. 


CONCLUSIONS 


To sum up, it seems justifiable to conclude that: 

1) Growing of Aspergillus niger on the surface of a foam film seems to be imprac- 
ticable. 

2) The method of spraying the inoculated medium into air may be feasible as a 
means of better aeration, but does not accelerate fermentation. 

3) The method of aeration and subsequent stirring with the aid of a vibro-mixer 
seems to be both convenient and practicable, reducing the duration of the process 
and achieving quite a high rate of conversion. Larger-scale experiments will 
be necessary, however, before this method can be properly evaluated. 
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LETTERS TO THE EDITOR 


A second order effect in plasticity. 


O. IsHAI AND M. REINER, Rheological Laboratory, Technion—Israel Institute of 
Technology, Haifa 


Williamson! found that when a cylinder of clay is rolled between two plates, a tube 
is formed. Egan and Jobling? repeated the experiment using different materials 
and found that the phenomenon is present with plasticine but is not present in 
polyisobutylene and silicone bouncing putty. We have used (1) cement with a water 
cement ratio of 26 percent (Figure 1), (2) cement mortar with a water cement ratio 
of 28 percent and a volume concentration of 42 percent of (standard Leighton Buz- 
zard) sand (Figure 2). In both cases tubes could be formed. In the second case 
pronouced radial cracks appeared. 

We suggest that this is a second order effect in the plastic deformation of solids. 
Polyisobutylene and silicone boucing putty are liquids, a solid being a material 
which has a yield point below which it does not flow. It seems that only granular 
materials are plastic solids. This would include polycrystalline metals where the 
single crystals form the granules. 

A second order theory will have to account for radial outward flow which mani- 
fests itself in both, the formation of a void at the centre and the appearance of circun- 
ferential tension causing radial cracks. 


REFERENCES 
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Figure | Figure 2 
Solid cylinder of cement paste shown on the left, Two tubes rolled from mortar paste. On the 
rolled into tubes as shown on the right. left view along axis, on the right—from side. 
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Analysis of thin rectangular plates, which are loaded by transverse concentrated 
forces 
Z. HAsHIN, Division of Mechanics, Technion—Israel Institute of Technology, Haifa. 


In the analysis of thin rectangular plates of elastic material, which are loaded by 
point forces, the loading is in general considered as a limiting case of either a surface 
load, distributed over a small rectangular or circular area, or a line load on a small 
interval. 

The solution is found in the form of infinite series. While the series, giving the 
deflection of the plate, converge quite rapidly, those by which the bending and 
twisting moments and shearing forces are expressed, converge very slowly. The 
convergence becomes ever worse as we approach the point of application of the 
load. Numerical evaluation of moments and shearing forces expressed by such 
series is all but impossible. 

This difficulty may be observed by considering the solution as being composed 
of two parts. The first part is the solution of an infinite strip, loaded by a point 
force, which may be expressed in closed form by a method due to Nadai! ; the 
second part is an infinite series which depends on the boundary conditions which 
have to be fulfilled in the specific problem under consideration. 

The bending moments and shearing forces may thus be given by a closed expression 
plus a very rapidly converging series, the convergence of which improves when 
approaching the point of application of the Joad. Numerical results may thus be 
obtained very easily for aJl points of a plate. 

The mathematical details will be published elsewhere. 
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